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Constraining the physical properties of Type II-P supernovae using 
nebular phase spectra 
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ABSTRACT 

We present a study of the nebular phase spectra of a sample of Type Il-Plateau supernovae with 
identified progenitors or restrictive limits. The evolution of line fluxes, shapes, and velocities 
are compared within the sample, and interpreted by the use of a spectral synthesis code. The 
small diversity within the dataset can be explained by strong mixing occurring during the 
explosion, and by recognising that most lines have significant contributions from primordial 
metals in the H envelope, which dominates the total ejecta mass in these type of objects. In 
particular, when using the [Oi] 6300, 6364 A doublet for estimating the core mass of the 
star, care has to be taken to account for emission from primordial O in the envelope. Finally, 
a correlation between the Ha line width and the mass of ^''Ni is presented, suggesting that 
higher energy explosions are associated with higher ^''Ni production. 
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1 INTRODUCTION 

Supernovae (SNe) are responsible for the production of most of the 
heavy elements that we observe in the Universe today. In a core- 
collapse SNe (CC-SN), this nucleosynthesis occurs both by hydro- 
static burning during the star's e volution and by explosive burning 
behind the shock after collapse. IWooslev & WeaveJ ( Il995i) anal- 
ysed the yields of elements that are synthesised in massive stars 
and concluded that stars of 15-25 M© are the typical stars responsi- 
ble for producing heavy elements. However, the initial mass func- 
tion shows that low mass stars are more common and therefore, 
depending on stellar yields, stars in the range 8-15 Mq could also 
make a non-negligible contribution to the chemical enrichment of 
the Universe. 

Type Il-Plateau (IIP) SNe are a subset of CC-SNe, which 
result from the core-collapse of massive H-rich stars, and there- 
fore display prominent H features in their spectra, along with 
an extended plateau of nearly constant luminosity in their light 
curves. They are the most comm on SN type m aking up ~40 per 
cent of all SNe per unit volume jSmartt et al.l l2009). Stellar evo- 
lutionary models predict th at stars in the range 8-30 Me should 



m 

explode as red supergia nts (iHeger et al."2003l: iLimongi & Chieffil 



l2003l : lEldridge & Tou^ |2004 ; Hirschi et al. 2Ci04h . However, the 
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progenitors of some IIP SNe have been identi fied, and apart from a 
few rare cases of blue supergiant progenitors jWalbometai]|l987l: 
Pastorello et al. 2005; Kleiseretal. 2011), they have been found 
to be red s upergiants in a lower than expected mass ra nge of 8-17 
Mo ( S martt et al.l2009l) . IUtrobin & Chugail ( l2008ll2009l) have mod- 
elled the light curves and spectral evolution of IIP SNe (some with 
progenitor identifications) and find larger masses of the exploding 
stars than those from the direct detections of the progenitors. This 
discrepancy between the progenitors obtained from direct imaging 
a nd those from mode llin g has also been analy sed and discussed 
in lSmartt et al.l ^20091) and lMaguire et al.l ( 1201 Oh using a sample of 
nearby IIP SNe. One way of studying this discrepancy is to inves- 
tigate and compare the nucleosynthetic yields from IIP SNe with 
identified progenitors. 

We do this by studying a sample of IIP SNe at epochs when 
we can see emission from newly synthesised heavy elements. Dur- 
ing the plateau phase of the light curve, which generally lasts for 
~ 80- 120 d, the light curve is powered byremission of the energy 
deposited by the shock wave, including H recombinations. At these 
times, the H envelope is optically thick and we can not see emis- 
sion from the metal core. However, after the H has recombined, the 
envelope becomes transparent and the core becomes visible. The 
light curve then enters the radioactive decay phase, which is pow- 
ered by the decay of ^*Co to ^^Fe. Radiative transfer effects are still 
significant in this phase, but decrease with time as the line opti- 
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cal depths fall. On the other hand, after approximately two years, 
the SNe generally become too faint to obtain good signal-to-noise 
spectra. Therefore, the optimal window for studying and interpret- 
ing the nucleosynthesis that has occurred is between one and two 
years post explosion. 

Very few IIP SNe have been studied in detail during the 
nebular phase and few late-time nebular spectra of IIP SNe ex- 
ist, with the notable exception of SN 1987A, which due to its 
proximity was very well monitored. This lack of late-time spec- 
tral data is partly due to the long exposure times that are re- 
quired to obtain sufficient signal-to-noise. The best observed IIP 
SNe at the necessary late-times ar e SN 1987A (e.g. , Phillips e t al. 
1990: Bouche t & panziger '1993: Me ikle et al.ll993l : [Wooden eXdX 
1993^) . SN I990E (iBenetti et al. 1994) . SN I999em 
200; 



.Leonard et al . 

Elmhamdi et alj l2003ah . SN 2004et jSahu et alj llOOd : 



Maguire e t al."2010'). SN 2007it ("Andrews etalj |201lh ."and SN 
2007od lAndrews et al. 2010; Inserra et al. 

It is generally assumed that the emission lines of specific ele- 
ments are emitted from regions containing predominantly that par- 
ticular element, i.e. that the O zone (formed in the late stages of 
pre-SN evolution or during the SN explosion) would be responsi- 
ble for the production of the O emission lines. Howev er, the final H 
envel ope mass for an 8-12 Mq star is ~6-8 M© (e.g. lHirschi et al.l 
|2004 . much larger than the mass of synthesised metals of ~0.5 Mq 
and it is not clear if emission by primordial metals in the envelope 
may dominate that o f synthesised ones. This has been shown by 
iLi & McCravl {l993§) to be the case in SN 1987 A for the [Can] 
7291, 7323 A and the Can 8600 A triplet, which are produced by a 
primordial abundance of Ca in the inner (velocity < 2500 km s"') 
region of the H and He envelope and not by newly synthesised Ca. 
Ca emission lines are very efficient cooling lines, and so even small 
quantities of primordial Ca can result in prom inent emission lines 
in the late-time spectra jpransson & Chevaliei|[i9 89). To obtain the 
observed flux of the Ca lines. iLi & McCravl 1 199331) needed only a 
Ca mass of 2x 10"'' Mo in ~5 Ms of H in the inner envelope of SN 
1987A, a factor of 10 times lower than that typically predicted for 
newly synthesised Ca in a SN 1987A-type progenitor star. 

Using mo dels that included realistic mixing of the nuclear 
burning zones, iKozma & FranssorJ h998ah showed that for SN 
1987A, the Can 8600 A triplet was at all times formed by emis- 
sion from primordial Ca in the envelope, whereas the [Can] 7291, 
7323 A lines were formed by both synthesised and primordial Ca 
until ~400 d, but mainly by primordial Ca after that. This does not 
imply that there is no newly synthesised Ca but instead that if it is 
present at late times, that the zones in which it is formed are too 
small and it cannot capture enough of the radioactive energy to cre- 
ate any strong emission, unless it is mixed with a grea ter mass of 
other elements such as H and He jLi & McCravlll993ah . 

Therefore, to calculate the relative emission strengths from 
primordial and synthesised material, the deposited energy in the 
different zones, and the dominant processes responsible for the pro- 
duction of these lines must be modelled using a spectral synthe- 
sis code. Then the synthesised mass of O, Ca, and Fe can be esti- 
mated by first subtracting the emission contribution from primor- 
dial material in the envelope. Previous studies of samples of CC-SN 
nebular phase spectra have mainly focused on stripped envelope 
SNe (Type lb (lacking H) and Type Ic (lacking H and He)) (e.g. 
Magdaetal. 20 08; Modjaz et al. 2008; Taubenberger et al. 2009; 
Milisavlievic et al.ll2O10h . and have studied the [Oi] 6300, 6364 A 



Ib/c SNe, focusing again on an analysis of the [O i] 6300, 6364 A 
doublet. Using the overall luminosity of the [O i] doublet and the 
ratio of the two components, they have estimated the synthesised O 
mass, assuming all the O luminosity is produced by freshly made 
O and not from primordial O in the envelope. We will in this paper 
show that this is a questionable assumption. 

In this paper, we study a sample of 25 late-time nebular phase 
spectra of nine IIP SNe with identified progenitors to investigate 
how their spectral features change as a function of time, and how 
the flux, profile shape and velocity of the features vary among SNe. 
This can then be linked to the material distribution and physical 
conditions in the SNe. We will apply a spectral synthesis model 
to the nebular phase spectra to study line flux ratios and interpret 
their evolution with time as well as their internal differences. This 
paper is complemented by Jerkstrand et al. (in prep.), which derives 
detailed nucleosynthesis results for one of the SNe in the sample, 
SN 2004et, by a complete analysis of the spectral evolution from 
the ultra-violet to far-infrared. 



profiles to look for possible signs of ejecta asp hericity in the line 
profile shapes and how they evolve with time. lElmhamdil i201lh 
have studied a sample of CC-SNe including IIP SNe as well Type 



2 SUPERNOVA SAMPLE OF NEBULAR SPECTRA 

As mentioned in Section [T] the current sample of IIP SNe with 
multi-epoch optical spectra at late times is limited to a small num- 
ber of objects. In particular we wanted to focus on SNe with in- 
formation available about their progenitor stars and ""^Ni mass es- 
timates. To increase this sample, observations were obtained of 
four IIP SNe at the 3.6-m New Technology Telescope (NTT) 
with the ESQ Faint Object Spectrograph and Camera version 2 
(EFOSC2) in low resolution spectroscopy mode and at the 8.2-m 
Very Large Telescope (VLT) with FOcal Reducer and low dis- 
persion Spectrograph 2 (FORS2). The SNe observed were SNe 
2007aa, 2008bk, 2008cn, and 2009N at epochs spanning -360- 
550 d post-explosion. In addition to the new data set, the literature 
was searched for late-time nebular spectra of IIP SNe with pro- 
genitor information (either a progenitor has been identified in pre- 
explosion images or an upper mass limit has been set) and data from 
the following SNe were included, SNe 1987 A, I999em, 2003gd, 
2004et, and 2005cs. The complete sample of spectra analysed here 
is detailed in Table[T] All the spectra have been corrected for Galac- 
tic and host galaxy extinction using the values given in Table|2]and 
they have also been corrected to the rest-frame of their respective 
host galaxies. The resolution of the spectra that are quoted were 
either taken from the literature or estimated from the full width at 
half maximum (FWHM) of the night sky lines. 

The ^*Ni masses and zero-age main sequence masses of the 
SN progenitors, where available, are given in Table|2] The proper- 
ties of the sample range from low-luminosity, low '*Ni mass events 
such as SNe 2005cs, 2008bk, and 2009N to higher-luminosity, 
higher ^^^Ni mass events such as SNe I987A, I999em, and 2004et. 
However, the range in main-sequence masses listed in Table|2]does 
not show as large a spread between over and under luminous events, 
with the masses clustering at the lower end of the range expected 
from stellar evolutionary theory (8-30 Mg). We will investigate 
here if there are any measurable differences between the spectral 
properties of the IIP SNe in the sample and then determine how 
these differences might relate to their progenitor and explosion 
properties. 
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Table 1. Spectral observations of the sample of IIP SNe studied, including epoch, wavelength range, resolution and source of the spectrum. The SNe for which 
new data are presented in this paper are noted as VLT+FORS2 and NTT+EFOSC2 as well as the grism with which they were taken, as discussed in Section|2] 



SN 


Epoch 


A RtingG 


Resolution 


Source 






(A) 


(A) 




1987A 


338 


3110-10930 


16 


Phillips et al. (1990) 




345 


3100-10610 


16 


PhilhDs et al. (1990) 




374 


5160-10690 


16 


Phillips et al. (1990") 




398 


3200-10220 


16 


Phillips et al. (1990") 




407 


3020-7780 


16 


Phillips et al. (1990) 




422 


3070-10740 


16 


Phillips et al. (1990) 




497 


3000-10990 


16 


Phillips et al. (1990) 




513 


3030-10970 


16 


Phillips et al. (1990) 




525 


3000-10990 


16 


Phillips et al. (1990) 




529 


4000-10350 


16 


Asiago Supernova Archive 




586 


3230-7500 


16 


Phillips et al. (1990) 




598 


3830-10170 


16 


Asiago Supernova Archive 


1999em 


312 


3350-10620 


21 


Elmhamdi et al. (2003a) 




319 


3410-7770 


21 


Elmhamdi et al. (2Q03a) 




391 


3500-9110 


10 


Elmhamdi et al. (2003^ 




465 


3360-10620 


20 


Elmhamdi et al. (2003a) 




510 


3560-9240 


10 


Elmhamdi et al. (2003a) 




642 


4200-9940 


28 


Elmhamdi et al. (2003a) 


2003gd 


493 


4200-9600 


13 


Hendry et al. (2005) 


2004et 


284 


4000-7820 


36 


Maauire et al. (2010) 




301 


4000-8450 


7 


Sahu et al. (2006) 




314 


4000-8450 


7 


Sahu et al. (2006) 




341 


3400-8020 


13 


Maguire et al. (2010) 




384 


3500-7750 


24 


Maguire et al. (2010) 




401 


4000-8450 


7 


Sahu et al. (2006) 




408 


4000-7730 


25 


Maguire et al. (2010) 




427 


4000-8450 


7 


Sahu et al. (2006) 




464 


4000-8450 


7 


Sahu et al. (2006) 


2005CS 


333 


5050-10470 


18 


Pastorello et al. (2009) 


2007aa 


376 


4300-9600 


10 


VLT-I-FORS2-I-300V 


2008bk 


524 


3500-9500 


14 


NTT+EFOSC2+gr#ll'' 




547 


3800-9500 


10 


VLT+FORS2+300V 


2008cn 


358 


3500-9500 


10 


VLT+FORS2+300V 


2009N 


366 


6000-10100 


14 


NTT+EFOSC2+gr#16'' 




406 


3500-9600 


10 


VLT+FORS2+300V 



"The resolution of the gr#16 and gr#14 gratings is very similar but gr#16 is centred at a longer wavelength than the gr#l 1. 

Table 2. Some of the derived properties for the sample of IIP SNe analysed in this paper including their extinction, distance, '*Ni mass and zero-age main 
sequence (ZAM) mass. 



SN 


E(B-V) 


Dist. 


Ni'* mass 


ZAM Source 






(Mpc) 


(Mo) 


(Mo) 



1987A 


0.19 


49.9x10^3 


0.075±0.005 


14-20 


1,2 


1999em 


0.06 


11.7±1.0 


0.042 


< 15 


3 


2003gd 


0.14 


9.3±1.8 


0.016 


7+6 


4 


2004et 


0.41 


5.9±0.4 


0.064 ±0.04 




5,6,7 


2005CS 


0.05 


8.4±1.0 


0.003 


8±2 


8, 9, 10 


2007aa 


0.014 


20.5±2.6 




< 12 


11, 12 


2008bk 


0.3 


3.9±0.4 


0.004 


Q+4 


13, 14 


2008cn 


0.35 


33.3±0.2 




13±2 


15, 16 


2009N 


0.15 


12.6±0.9 


0.009±0.004 


<16 


16, 17 



R EFERE NC ES - (l)LP hillips et alj Jl990l): (2) Smartt 1 2009); (3) Elmhamdi et al. 1 2003a) (4) Hendrv et al. 1 2005); (5)Sahu et al.' ilOOf^: (6)'Magui re etal] 

ilOldt): (7)ICrockett et al. t201ll); (8)lMaund. Smartt & Danziger (2005) ( 9) .Li et al.(.2006)) (10) P asto rello et al. (2009); (1 D Smartt et al. I 2009J ; (12) 
IChomock et aljj2010l) ; (13) lMattila et al.H2008l) ; (14) lvan Dvk et alj j2010l) ; (15) lElias-Rosa et a\WW^ : (16) Fraser et al. (in prep.); (17) Takats et al. (in 

prep.) 
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3 DESCRIPTION OF RADIATIVE TRANSFER MODEL 

Full details of the spectral synthe sis model used to study the nebu- 
lar phase spectra can be found in ljerkstrand et alj ( 1201 ll) and Jerk- 
strand et al. (in prep.). Here we describe only the basic setup of the 
model and its important input parameters. The model uses as an 
input the 1 2 Mp explosion model o f a solar metallicity progenitor 
taken from lWoosley & Heged ( 120071) . A 12 Mq model is chosen be- 
cause this mass is typical of the progenitor mass estimates obtained 
from pre-explosion images of the sites of IIP SNe and stellar evo- 
lutionary models (Smartt et al. 2009, and references therein). The 
late stellar evolution and explosion of stars less massive than 12 
Mq involve significant uncertainties, and few if any useful models 
are available. In Jerkstrand et al. (in prep.) we discuss models for a 
wider set of masses. 

The ejecta is divided into a spherically symmetric core and a 
spherically symmetric envelope, with the core containing the zones 
of the heavier elements (Fe/He, Si/S, 0/Si/S, 0/Ne/Mg, 0/C), along 
with 15 per cent of the mass of H and 60 per cent of the mass of 
He mixed in from the envelope. The zones making up the core re- 
gion are broken into clumps (100 for each zone), with the clumps of 
each zone having the same density and being randomly distributed 
throughout the core region. This simulates macroscopic mixing in 
the core and removes the rigid shell-like structure that is known 
not to be an accurate representation of the ejecta structure post- 
explosion, as will be discussed in Section[4l We allocate t he filling 
factors using the same method as in ljerkstrand et al.l ( 1201 ll) . assum- 
ingthe same density for the synthesised O as derived for SN 1987A 
byU^i & McCrav ( 1992). At epochs in the range 100-1000 d post- 
explosion, the main source of luminosity is due to the radioactive 
decay of '^Co and its mass can be constrained by measuring the 
bolometric luminosity during the radioactive tail phase of the SN. 
For this 12 Mq model, a '*Ni mass of 0.064 Mg and a core veloc- 
ity of 1800 kms"' were chosen to agree with the values obtained 
from observations of SN 2004et. The model computes the deposi- 
tion of gamma rays and positrons, followed by iterative solutions 
to the equations of st atistical and thermal eq uilibrium, as well as 
the radiation field. See lJerkstrand et ^ j201 Ihfor more details. We 
note that for low velocity/low '^Ni mass SNe such as SN 2008bk, 
the conditions for the steady state calculation assumed in the model 
may not be met at all times, and the model line fluxes are therefore 
of limited accuracy for these types of SNe. 



4 LINE PROFILE ANALYSIS 

Standard stellar evolution models of massive stars (e.g. 
Wooslev & Weaver 1 19951 : iThielemann." Nomoto & Hashimoto! 
19961 : iHirschi et alj2004h describe how the pre-explosion structure 
of the progenitor has an ordering of layers of elements from each 
nuclear burning stage, with the outermost layer being H and the 
innermost Fe. In the absence of mixing, this onion-layer structure 
would hold for SNe after explosion and so if the emission of 
each particular element is formed in a different shell containing 
predominantly that element, then there would be a correlation 
between the elemental mass numbers and the width of their 
emission lines. The lightest elements (H, He) would have broad 
lines with pronounced flat-tops, and the heaviest elements (Si, Si, 
Fe, Ni, Co) would have narrow lines. At the other extreme, if the 
elements were distributed uniformly throughout a sphere, then the 
line profiles would be parabolic in shape. 

Studies of the observational properties of SN 1987 A have 



found very similar line widths for a wide range of elements 
with no flat-tops, which suggests extensive mixing of the ejecta 
dFransson & ChevaUeJ 1 19871 : iMeikle et alj 1 19891 . Il993l) . This is 
supported b y the early detection of gamma rays an d X-rays from 
SN 1987A ('Po tani et alJll987l : ISunvaev et alJI 19871) . which would 
only occur if ^*Co and ^^C o were mixed out i nto th e inner envelope 
in co ncentrated clumps dPinto & WooslevI Il988l : iKumagai et al.l 
Il989l) . The detecti on of a hig h excitation line. He 1 10830 A at times 
greater than 10 d telias et a l. 1988) could also only occur if it was 
reionised by gamma r ays emitted in the He zones from outward 
mixed, decaying ^"^Co jGrahadll988l : lFassia et aljl998l) . 

4.1 Line profile shapes 

Here we investigate the shapes of the profiles of the nebular phase 
emission lines in multi-epoch spectra for our sample of SNe, to 
study the distribution of elements in velocity space, and determine 
if signatures of mixing are also seen in the spectra of our sam- 
ple. Figure[T]shows the velocity profiles of the prominent emission 
lines present in SN 2004et at 300-464 d, while Figure|2| shows the 
velocity profiles for SN 1999em in the range 312-510 d, and Fig- 
ure|3]shows the profiles for low-luminosity SNe 2008bk (524, 547 
d) and 2009N (366, 406 d). The lines shown are Ha, [Oi] 6300, 
6364 A, [Can] 7291, 7323 A, and [Feii] 7155, 7172 A. Each of the 
lines shown has had a linear continuum subtracted and has been 
normalised to its peak value. The [Oi], [Can], and [Feii] lines are 
doublets and to compare the emission line shapes, the profiles have 
been centred on the rest wavelength of the blue components of the 
doublets. 

None of the emission lines of these SNe shown in Figures [Tj 
[2] and [3] has a flat- topped profile at any epoch and the lines are in- 
stead roughly Gaussian in shape. However, the convolution of the 
emission lines with the point spread function (PSF) of the instru- 
ment could result in a roughly Gaussian shape being seen for an 
intrinsically flat-topped profile if the intrinsic width of the line is 
smaller or comparable to the resolution of the instrument. To in- 
vestigate the effect, we have convolved a flat-topped profile (with 
the width of the Ha line) with a Gaussian (with the width of the 
instrumental resolution). However, before doing this, the observed 
FWHM (FWHMobs) must be corrected to account for the resolution 
of the instrument (FWHM|ns,i.) as measured from the FWHM of the 
night sky lines in each of the spectra to obtain the original width 
of the line before convolution with the telescope instrument. The 
corrected FWHM is given by, 

FWHM,,,, = ^FWHMl^^ - FWHMl^^^ (1) 

Then for each SN in the sample, we convolved a flat-topped line 
profile (with the width of the Ha line, corrected for the instrumen- 
tal resolution) with a Gaussian (with the width of the instrumental 
resolution) to determine if any of the line profiles were flat-topped. 

To demonstrate this, in Figure U a comparison between the 
convolution of the resolution of the instrument with both a flat- 
topped and Gaussian profile is shown for two different FWHMcoi,- 
values. In the top panel of Figure |4] using a FWHMcon of 20 A, 
it is not possible to distinguish from the output between the two 
scenarios. However, in the bottom panel with a FWHMcoir value 
of 28 A, the result is different, with the flat-topped profile giving a 
broader, flatter-topped line profile. The FWHM at which it becomes 
difficult to distinguish between the intrinsic profile types for the 
quoted resolutions is a FWHM of -29 A for the NTT+EF0SC2 
(resolution of 14 A) and a FWHM of -22 A for the VLT-I-FORS2 
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2004et 40' d 




-7500 -3750 3750 7500 

Velocity (km s"') 



2004et 408 d 




Velocity (km s"') 



2004et 42? d 




2004et 464 d 




Figure 1. Line profiles for SN 2004et between 300 and 464 d. The [Oi], [Can] and [Fen] lines shown are doublets and have been centred on the blue 
component, as discussed in Section l4~n The profile of Ha is shown in black, [Oi] 6300, 6364 A in red, [Can] 7291, 7323 A in green, and the [Fen] 7155, 
7172 A doublet is shown in blue. 



(resolution of 10 A) at the position of Ha assuming a signal-to- 
noise of greater than ~20. 

We have found for all the SN spectra in our sample that it 
is possible to distinguish between these two scenarios and so the 
roughly Gaussian shaped profiles that are seen are found to be in- 
trinsic to the SN emission, and this immediately suggests that mix- 
ing of the elements must occur in the SN ejecta to allow elements 
to be at zero velocity. Only for the spectrum of the narrow-lined SN 
2008bk obtained at the NTT+EF0SC2 with a resolution of 14 A, 
was it found not to be possible to distinguish between the intrinsic 
Gaussian profile and a flat-topped one. However, another spectrum 
of SN 2008bk was obtained at the VLTH-FORS2 with a resolution 
of 10 A, where it was possible to distinguish between these two 
scenarios and the line profile of Ho- of SN 2008bk was found not to 
be intrinsically flat-topped. 

In Figures[T]|2]and[3] it should be noted that the profiles of the 
different lines for a particular SN appear to have similar line profile 
shapes, which again suggests that the different elements produc- 
ing the emission lines have a similar spatial distribution in the SN 
ejecta. In the early spectra of SN 2004et in Figure [T] and for SN 
1999em in Fig.|2l it appears as if the blue wing of the [O i] 6300 A 
line extends to higher velocities than the other lines. However, this 
early blue wing of the [Oi] 6300 A line disappears at later times 
and is most l ikely caused by a blend w ith another emission line. 
In particular, lOessart, Livne & Waldmani ( 120 lOi) showed that there 
are many weak Fe lines i n this region that could b e blending with 
the [O i] doublet (see also lFransson & Kozmall2002[) . If the H emis- 
sion lines were being produced solely in the envelope and the other 
elements from a central region then it would be expected that the 
wings of the H profiles should extend to greater velocities, which 
does not appear to be the case. However, the velocities that can be 
measured from these emission lines are wavelength dependent and 
also must be corrected for the resolution of the instrument used to 



obtain the observations and therefore, a more quantitative analysis 
of this is presented in Section l4~2l 

An evolution with time in the line profile shapes could be ex- 
pected if at early times the emission mainly comes from regions 
close to the ^''Ni but then as the gamma rays propagate outward, 
emission from other regions is seen. However, no evolution is seen 
with time, which could suggest that ^^Ni is not concentrated in the 
core of the SN but instead is distributed by mixing to regions farther 
out in the ejecta. 

Asymmetry of some of the line profiles is seen in Figures[T]|2] 
and|3] The nebular phase [Or] 6300, 63 64 A doublet of SN 2004et 
was investigated bv iutrobin & Chugail ( 12009 ) for signs of possible 
asymmetry in the ejecta. They fitted three Gaussian components to 
each of the lines of the [O i] doublet and inferred a bipolar structure 
of the line-emitting gas in the inner layers of the SN envelope. They 
noted that this asymmetry refers to the line-emitting gas, which is 
not identical to the overall O distribution. An attenuation of the red 
wing (blueshift) of the line pro files of Hg and [Oi] i s seen for SN 
2004et at -300 d and shown in lMaguire et al.l ilOld) to be coeval 
with other indicators of dust formation. Dust condensa tion within 
the m etal-rich ejecta of SN 1987A was first noted bv ILucv et al.l 
( Il989t) from the asymmetry of optical emission lines. Signatures 
of dust f ormation such as line a ttenuation were also seen for SN 
1999em dElmhamdi et alj2003al) at -500 d post explosion, and SN 
2009N appears also to show a blueshifted profile for the [O i] 6300 
A line at both the epochs shown in Figure |3] which is most likely 
due to dust formation in the ejecta. 

Other signs of dust formation include a decrease in lumi- 
nosity at optical wavelengths accompanied by an increase at NIR 
wavelengths. Signatures of dust formation have been seen at mid- 
infrared wavelengths of IIP SNe to occur at similar epochs to the 
appea rance of t he blueshifted profiles for SNe 1987 A dMeikle et al.l 
Il993h . 2003gd dSugerman et al.ll2006l : lMeikle et alj|2007l) . 2004di 
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Figure 2. Line profiles for SN 1999em between 312-510 d. The lines are 
as specified in Figure^ 
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Figure 3. Line profiles for two low luminosity SNe 2008bk and 2009N. The 
lines are as specified in Figure[T] 



( lMeikleetal.l201lh . and 2004et jKotak et al.l2009h . Therefore, the 
shapes of the emission lines along with other observational param- 
eters, tell us not only about the velocity distribution of the ejecta but 
also about the onset of dust formation and any possible asymmetric 
distribution of elements in the ejecta. 



4.2 Velocity distribution 

To study the distribution of elements in the ejecta and velocity at 
which the bulk of the emission occurs, the FWHM of the promi- 
nent emission lines in the sample of SNe were measured by fitting 
Gaussian profiles to each of the lines. The FWHM used were cor- 
rected as discussed in Section HTl for the resolution of the instru- 



Figure 4. Comparison of the convolution of the instrumental resolution ob- 
tained with the VLT+FORS2 of 10 A (measured from night sky lines) with 
both a flat-topped and a gaussian profile. In the top panel, a FWHMcon of 
20 A is used, which corresponds to a post-convolution value of ~22 A and it 
can be seen that the output line profiles are very similar and it is not possible 
to distinguish between them. In the bottom panel, an input FWHMcon of 28 
A is used, which corresponds to a post-convolution value of ~30 A and the 
output fine profiles for the flat-topped and gaussian inputs are found to be 
different. 



ment used to obtain each spectra using Equation[T] From these cor- 
rected FWHM values, the half width at half maximum (HWHM) 
velocities can be obtained to determine the velocity to which the 
different elements are found in the ejecta with respect to the zero 
velocity of the line profile. 

The HWHM velocities of the sample of SNe studied here are 
shown in Figure[5] The SN sample is split into three groups based 
on the average velocities of the emission lines. For the [Oi] and 
[Ca ii] doublets, the average value of the measurements of the two 
lines was used. Limits are shown for some of the emission lines 
when the individual components were not well enough resolved to 
fit Gaussians to each of the lines and were obtained by fitting a sin- 
gle Gaussian to the whole line profile. Therefore, these limits are 
very conservative upper limits for the widths of the individual lines. 
The outer wings of the line profiles could also be calculated from 
either the position at which the line flux goes to zero or the loca- 
tion of the first local minimum outward from the centre. The posi- 
tion of these outer wings would then tell us the velocity to which 
a certain element extends. However, the uncertainties in measur- 
ing these values are large due to the effects of blends of lines, 
contributions from the continuum and noisy spectra where the lo- 
cal minimum does not represent the tru e line minimum. Therefore, 
like iDessart, Livne & Waldmanl i2010l) . we have decided to use the 
more reliable HWHM for calculating the velocities, since it is more 
easily measured and compared among SNe. 

Extensive studies of the emission line profiles of SN 1987 A 
have been performed in the past and as a consistency check, we 
have compa red our measuremen ts of SN 1987A to those of previ- 
ous studies. lMeikleetal.l ( ll993h studied 1-4 |j.m spectroscopy of 
SN 1987A and found that the Oi 1.129 iJ.m velocities were lower 
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than those of H, in agreement with Figure |5] The [Fe ii] and [Ca ii] 
velocities of SN 1987 A have only upper limits plotted, but the 
[Can] lines are consistent with having been produced i n the in - 
ner regions of the H envelope as detailed in, Li & McC rayj l ll993ah . 
The [Oi] velocities a re also consistent with t h ose sum marised in 
iLi & McCravl ( Il992h . iLi. McCrav & SunvaevI h993b l) found that 
the Fell, Con, and Nin Unes of SN 1987A have FWHM of 3000- 
3500 km s"' , which are consistent with the HWHM velocity limits 
of [Fell] shown here. The measurements of SN 1987A are consis- 
tent with those of previous studies, which suggests that our method 
of fitting the Gaussian profiles is sound. 

In Figure |5] there is a much larger difference between the Ha 
and [O i] velocity for SN 1987A than for the other SNe in the sam- 
ple. This could suggest that the H envelope plays a role in con- 
tributing to the flux seen in the prominent emission lines for the 
other SNe, since the velocities for the different elements are rela- 
tively similar. The immediate physical interpretation of this is that 
a major fraction of the [O i] emission does not come from newly 
synthesised O, but from primordial O that is present in the H zones 
or from O clumps mixed to high velocities. For the low-luminosity, 
low velocity events (SNe 2005cs, 2008bk, and 2009N) shown in 
Figure [5] actual measurements (instead of limits) of the expansion 
velocities were obtained for all the lines and it can be seen that Ha 
has higher velocities than the other lines, although the difference is 
not as great as in SN 1987 A. 

In iDessart. Livne &" WaldmanI jlOlOt) , radiation- 
hydrodynamic simulations of IIP SNe were performed. They 
used a ID model without mixing and suggested that their models 
of the chemical distribution of the ejecta in velocity space can 
be used to place strong constraints on the main-sequence masses 
of IIP SN progenitors. As shown above, the line profile shapes 
indicate that mixing is likely to have occurred in the ejecta of IIP 
SNe, and this argues against t he assumption of chemical s tratifi - 
cation in the SNe. While the iDessart, Livne & WaldmanI ( 1201 Oh 
models are an informative first step in modelling nebular spectra 
to probe progenitor and explosion characteristics, the mixing of 
material in velocit y space is likely an import a nt par ameter for 
future exploration. IDessart. Livne & Waldman ( 120101) used the 
width of the [O i] 6300, 6364 A doublet to estimate the velocity 
of the 0-rich zones as it w as shown for SN 1987A at ~500 d by 
iKozma & FranssorJ l [l9983) that the [O i] lines were formed in the 
O-rich zones. However, in Section lSTTI we describe how a radiative 
transfer model can be used to determine the contribution to line 
emissions from different zones in the SN ejecta and it is shown 
for a 12 Mo model that the flux in the [Oi] 6300, 6364 A doublet 
has approximately equal contributions from primordial O in the 
H-rich zones and synthesised O in the O-rich zones. However, as 
the mass of the progenitor increases, as in the case of SN 1987 A, 
the contribution from the O core becomes dominant. 

To investigate the effect of O emission from the envelope on 
the measured FWHM of the [Oi] 6300, 6364 A doublet, we used 
the spectrum resulting from the model detailed in Section [3] We 
compared the FWHM values obtained from the model spectrum 
with and without the contribution from the H zones and found that 
by excluding the contribution from the H zones, the FWHM of the 
[Oi] 6300, 6364 A doublet in the model spectrum was reduced by 
~8 A, which corresponds to a HWHM velocity difference of -200 
km s"' . This would suggest that using the HWHM of the [O I] 6300, 
6364 A lines as an indicator of the extent of the spatial distribution 
of O would result in small but measurable overestimates of the ve- 
locities. The relatively small reduction of ~8 A when excluding 
the contribution from the H zones is due to the fact that while half 



Table 3. Primordial and .syntliesi.sed masses of tlie listed elements from the 
12 Mq model used in this paper, along with the ratio of primordial to syn- 
thesised mass for each of the elements. The Fe mass is calculated at 400 d 
for a "^'Ni mass of 0.064 Mq. 



Element 


Synthesised 


Primordial 


Synthesised/Primordial 




Mass (Mq) 


Mass (Mq) 


Mass Ratio 


O 


3.08E-01 


4.45E-02 


6.9 


Si 


4.12E-02 


7.12E-03 


5.8 


S 


2.63E-02 


3.65E-03 


7.2 


Ca 


2.40E-03 


6.41E-04 


3.7 


Fe 


6.5E-02 


1.22E-02 


5.4 



of the [O i] flux is emitted by these zones, the main H zones con- 
tributing to the [O i] flux are those that are mixed into the core and 
therefore, have a lower velocity, than the average of all the H zones. 



5 LINE FLUX ANALYSIS 
5.1 Model results 

The synthesised and primordial masses from our 12 Mq model are 
given in Table[3]for some of the main elements, along with the ra- 
tio of primordial to synthesised mass for each of these elements. It 
can be seen that the synthesised masses are all larger than the pri- 
mordial masses. However, when it comes to the observed spectrum 
it must be remembered that the flux in the emission lines will also 
depend on how efficient an emitter the element is per unit mass and 
how the energy is deposited in the ejecta. 

The spectral synthesis model can be used to calculate the con- 
tributions to the luminosity of the prominent emission lines ob- 
served during the nebular phase from the various regions of the 
SN ejecta. In this section, we compare the output of the model to 
the observations in the range of 300-650 d. We make specific in- 
depth comparisons at 400 d, when most of the SNe in the sample 
have observed spectra. At this epoch, our model with a ^*Ni mass of 
0.064 Mq, gives a total deposited energy of 1.9 x lO"*" ergs, which 
is ~75 per cent of the emitted energy (-25 per cent of the gamma 
rays escape). The main characteristics of the model at 400 d are 
detailed in Table O including the energy deposited in each of the 
zones as well as the channels through which this deposited energy 
is released (heating of the gas, ionisations and excitations) and the 
main cooling lines. In this table, F, is the fraction of the energy de- 
posited into zone i, fj,^^, fl^^^ and are the fractions of this energy 
going into non-thermal excitations, ionizations and heating, respec- 
tively. Finally, /j is the fraction of the cooling done by a given line 
in a given zone. 

It can be seen that approximately 2/3 of the total energy is 
deposited in the H zones, with the next highest amount of deposited 
energy (-1/6) going into the He zone. Since most of the energy is 
deposited in the H zones and mainly released by heating of the gas, 
it is expected that nebular phase spectra should be dominated by 
thermal line emission from H-rich material. The main cooling lines 
in the H zones during the first year are lines at UV wavelengths such 
as the Mgii 2795, 2802 A, and Feii 2382, 2600 A lines. However, 
there is also an appreciable contribution from cooling lines in the 
optical that can be analysed. Although the He zone gets -1/6 of 
the energy, it does not produce any distinct optical lines that can be 
studied. For most of the zones at 400 d, it can be seen that heating 
is the dominant channel for the deposited energy, with ionisation 
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Figure 5. The HWHM (half the FWHM) velocity of the emission hnes of Ha (black), [Oi] (red), [Feii] (blue), and [Can] (green) of SNe 1987A and 2008cn 
(top panel), SNe 1999em, 2004et, and 2007aa (middle panel), and SNe 2008bk, 2009N, 2005cs, and 2003gd (bottom panel) are shown. The SNe are grouped 
based on the average velocities of the emission lines. These velocities have been corrected for the intrinsic resolution of the instrument used to obtain each 
spectrum. 



the next most important, while only a small percentage of energy 
goes into excitations. 

The excitation of the lines can in general be divided into those 
excited thermally and those excited by recombinations. The recom- 
bination lines arise as a result of the non-thermal ionizations and in 
some cases, like H, photoionisations. They are therefore insensi- 
tive to the temperature, and mainly depend on the efficiencies for 
ionisation by the gamma-rays and the ionising UV radiation. The 
thermally excited lines are responsible for most of the cooling of 



the ejecta. The strength of a given line from a specific abundance 
zone is determined by the fraction of energy deposited in that zone 
and the relative abundances in this zone. However, the efficiencies 
of different elements in cooling are very diflrerent, depending on 
their atomic structure. As an example the [Can] 7291, 7323 A lines 
are typically a factor of 10^ - 10^ times more efficient per atom 
in cooling than the [Oi] 6300, 6364 A lines. A small abundance 
of Ca can therefore dominate that of the much more abundant O 
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Figure 6. The contribution to the Ha luminosity from non-thermal ionisa- 
tions, /|^„"' (solid line) and from photoionisations, fl^f'' (dashed line). 



jpransson & ChevalieJI 1989b . The luminosity, L(/i) of the different 
lines can be estimated from Table |4]using. 



L(A) = E,„, F, fl,^^,, f- 



(2) 



where B,o, is the total energy deposition (in the whole ejecta) and 
the other quantities are defined above and in Table|4] We discuss the 
individual lines here, referring specifically to the model described 
above. 



5.1.1 Ha 

Ha is not an efficient cooler because the n=3 state has too high an 
excitation energy (10.2 eV) and Ho- is instead created by recom- 
binations, and to a lesser extent by non-thermal excitations. The 
importance of phot oionisations from the n= 2 state by UV photons 
was first shown by Kirshner & Kwanl l ll975h . The model results for 
the relative contributions are sho wn in Figure [6| a nd found to be in 
good agreement with the results of lXu et al] ( ll992h . is the frac- 
tion of the energy deposited in the H zones going into non-thermal 
ionisations and '' is the fraction of the energy deposited in the 
H zones going into Balmer photoionisations. The value 0.14 is the 
fraction of non- thermal ionisation energy reemitted as Ha (1.89 
eV/13.6 eV) and 0.56 is the fraction of Balmer photoionisation en- 
ergy reemitted as Ha (1.89 eV/3.4 eV), where 1.89 eV is the energy 
of Ha photon and 3.4 eV is the ionisation potential of the Balmer 
continuum. Using our calculation, which includes a detai led com- 
putatio n of the internal UV field , we co nfirm the results in lXu et al] 
( 1 19921) and iKozma & FranssorJ ( Il992h . that Balmer photoionisa- 
tions are dominant in the H ionisation rates up to ~400-500 d. At 
this time a transition to the non-thermal dominated phase occurs, 
when the Balmer continuum becomes optically thin and direct ion- 
isations from non-thermal electrons dominate. At 400 d, we find a 
deposition- weighted ratio of Balmer to non- thermal ionisation of 
~3 (see Figure|6j. 



5.1.2 [Oi] 6300, 6364 A 

The emission in the [O i] 6300, 6364 A lines is caused by a com- 
bination of emission from synthesised and primordial O. There are 
three zones that contain O as their main component; the 0/Si/S, 
0/Ne/Mg, and the 0/C zones. In the 0/Si/S and O/C zones, mod- 
elling has shown that molecules of SiO and CO form and ar e 
responsible for nearly all of the cooling l lLiu&Dalgmiolll995[) . 
This is supported by the detection of SiO and CO in SN 1987 A 
dSpyromilio etalj|l98 8': Aitken et al. 1988), as w ell as in other II P 
SNe such as SN 2004et (Kotak et al. 2009; Ma guire et ai]|2010l) . 
Therefore, we can expect to see thermal emission from synthesised 
O from the 0/Ne/Mg zone only. In this zone, the [O i] 6300, 6364 
A doublet is one of the dominant cooling lines doing about half of 
the total cooling and provides a significant contribution to the [O i] 
6300, 6364 A lines luminosity. 

However, there is also a contribution to the [O i] doublet from 
primordial O in the H zones, where the doublet is found to do a 
few per cent of the total cooling. From Table |4] the fraction of the 
total emission in the [O i] lines from the 0/Ne/Mg zone is ~ 0.025 x 
0.52x0.55 ~ 0.0072, while it is ~ 0.67x0.35x0.044 ^ 0.010 in the 
H envelope. The low fraction of the cooling done by the [O i] lines 
in the H envelope is therefore balanced by the large total energy 
deposition in this zone. 

This leads to the important conclusion that in IIP SNe with 
progenitor masses of ~12 M©, the [Oi] emission is formed by ap- 
proximately equal contributions from primordial and synthesised 
O. This is perhaps surprising since there is ~6 times more synthe- 
sised than primordial O present (assuming solar metallicity), but 
this is negated by the fact that the H zones absorb ~8 times more 
gamma ray energy than the O zones, and also that the O/Si/S and 
O/C zones are cooled by molecules. 

This is different to the results that were obtained for SN 
1987 A using an 20 M© model, where it was found that [Oi] 
e mission lines are formed i n the zones containing syn thesised 
O jKozma & Franssonlll998al) . IWooslev & Weaved jl995h showed 
that the O mass depends sensitively on the progenitor mass; a 20 
Mo main-sequence star would produced an O mass of 1.8 Mq and 
the synthesised O mass would be ~40 times greater than the pri- 
mordial O mass, and so the contribution of synthesised O to the 
total [O i] 6300, 6364 A lines should dominate. These models for 
SN 1987 A also assumed a lower amount of O in the envelope due 
to CNO mixing than obtained in the stellar evolution models used 
here. 



5.7..? [Can] 7291, 7323 A 

Synthesised Ca exists mainly in the Si/S zone (2.1 x 10"^^ Mq), but 
also in the Fe/He zone (3.2 x 10""* Mq), while the total mass of pri- 
mordial Ca present in the envelope is 6.4 x lO""* Mq( ~25 per cent of 
the synthesised mass). As was previously shown by ' Li & McCravl 
( 1993a) and Kozma & Fransson ( 1998a) for SN 1987A, it is found 
for our 1 2 Mq model that the [Ca ii] 729 1 , 7323 A doublet is formed 
nearly completely by primordial Ca in the inner regions of the H en- 
velope. The [Can] 7291, 7323 A emission from the Si/S and Fe/He 
zones is no more than ~ 10 per cent of the total flux at any time. This 
is because the masses of these zones are small compared to the H 
zones and so receive only a small amount of gamma-ray energy. 
At later times, the [Ca ii] lines will receiv e important contributions 
from fluorescence ( lLi&McCravlll993ah but at -400 d the [Can] 
emission from the H envelope is dominantly thermal. 



© 0000 RAS, MNRAS 000, 000-000 



10 K. Maguire et al. 

Table 4. Output values for the model at 400 d: total energy deposition, energy channels, and dominant cooling lines for each zone. 



Zone 


Fe/He 


Si/S 


O/Si/S* 


O/Ne/Mg 


O/C* 


He" 


W 


Energy deposition % (.^-one^ 


7.5 


1.9 


2.5 


2.5 


2.9 


16.2 


67 


Tir^>nlfrnn\MT} /'if I^n^vov r^nnytnflv 
















Heating %" {f~"/J) 


76 


52 


43 


52 


43 


43 


35 


lonisation % (ij^^) 


18 


29 


47 


42 


49 


44 


37 


Excitation % (/j;"^"'') 


6 


19 


10 


6 


8 


13 


28 


Breakdown of Cooling Lines: 
















[Can] (7291, 7323 A) 


2.4 


47 





7.2 





2.3 


27 


[0 1] (6300, 6364 A) (/g^"™) 











55 








4.4 


[Fen] (7155,7172A)(/,'°5"p 


7.4 


0.6 





0.3 





0.3 


1.7 



Energy deposition channels are calculated for the H and He zones using a mean, weighted by the energy deposition in each zone. 
''The zero values are due to molecule formation, which we assume is responsible for the cooling in these zones. 
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Figure 7. From top to bottom, the ratio of the flux in the [Ca 11] doublet to Ha and the ratio of the [O i] doublet to Ha. The dashed black line shows the output 
from the model. 



5.1.4 [Fen] 7155, 7172 A 

The [Fen] 7155, 7172 A doublet has a flux ratio of its components 
of ~3 to 1 if the parent multiplet is in local thermodynamic equilib- 
rium and optically thin. The state from which the doublet comes is 
only 1.95 eV above the ground state, so this transition is thermally 
driven. The [Feii] 7155, 7172 A doublet is formed by contributions 
from both synthesised Fe from the Fe/He zone and primordial Fe 
in the H envelope. From Table|4] one finds from the products of the 
energy deposition, heating fraction and cooling fraction, that there 
are equal contributions from the Fe/He zone and the H envelope to 
this line. This line, like the [O i] and [Ca n] doublets, is therefore a 



mix of synthesised and pri mordial Fe, vyhich is agreem ent with the 
results for SN 1987A (e.g. lKozma & Franssonlll998ah . 

5.2 Application to observed line flux ratios 

Having obtained a theoretical interpretation of the origin of the flux 
in the prominent emission lines and the dominant processes in- 
volved (ionisation in the case of Ha and heating for the rest of the 
studied lines), these results can then be compared to the values ob- 
tained from our late-time nebular phase spectra using flux ratios of 
the prominent emission lines. Line ratios are a convenient method 
of comparing emission line properties of SNe both as a function of 
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Figure 8. Ratio of flux in [Ca ii] to [O i] and ratio of flux in [Fe n] to Ha. Tlie dashed black line sliows the output from the model. 



time and also between different objects, because any uncertainty in 
calibrating the absolute flux of the spectra is removed. 

The line ratios were measured by fitting a Gaussian to the 
line profiles and then integrating the flux under the curve and 
subtracting ofi' a linear continuum. In the case of the doublets of 
[Oi], [Can], and [Feii] where the individual component lines were 
blended, both lines were fitted simultaneously with two Gaussians 
to obtain the total flux. In Figures[7l[8]and|9l the ratios of the promi- 
nent emission lines (Ha, [Oi] 6300, 6364 A, [Can] 7291, 7323 A, 
and [Fen] 7155, 7172 A) for the nine SNe described in Section|2] 
are shown, along with the model curve obtained from the spectral 
synthesis model detailed above. The flux ratios outputed from the 
model are for a 12 Mq model with a core velocity of ~1800 kms"'. 
As noted in Section [3] this model is not suitable for calculating 
line-flux ratios of low-luminosity, low-velocity events such as SN 
2008bk (core velocity of ~600 km s"' ) due to strong adiabatic cool- 
ing at these velocities, which makes steady-state temperature cal- 
culations unreliable. Therefore, comparisons of low- velocity events 
such as SNe 2008bk and 2009N with the model output in Figures 
|7l[8]and|9]should be treated with caution. 

Firstly, it can be seen in Figures |7] and [8] that there is a trend 
of increasing flux from the [Can], [Oi], and [Fen] lines with time 
compared to Ho-. This can be understood based on the transition, 
as detailed in Section [5. 1.1 1 of the Balmer continuum to an opti- 
cally thin regime and shown in Figure |6] The contribution of the 
B aimer-continuum photoionisation to the total ionisation rate de- 
creases with time and hence, the Ha luminosity also decreases. 
Also, we see in Figure [10] that the fraction of the total luminosity 
that is carried by the optical cooling lines increases up to 500-600 



d, as the temperature falls to favour optical over UV line cooling. 
At ~600 d, infrared lines take over and the optical lines decrease. 

The [Can] 7291, 7323 A to Ha model ratio as a function of 
time is plotted in Figure [7] and found to be slightly lower at early 
times, but is in reasonable agreement with the observed flux ratios. 
The [Ca n] to Ha ratio displays little variation between the different 
SNe and this can be understood by both of the lines being formed 
in the H zones, and therefore is independent of its mass, and the 
relative amounts of energy going into heating (responsible for the 
[Can] emission), and ionisation (responsible for the Ha emission) 
are not very sensitive to the density and temperature. 

The [Oi] 6300, 6364 A to Ha model ratio shows the same 
trend as the observed flux but is lower at all epochs. Since this is the 
same underproduction trend as seen at earlier times for the [Ca n] 
to Ha ratio. Ha may be overproduced in the model. The scatter 
in the ratio between different SNe does not seem to show any trend 
with '*Ni or progenitor mass. Although the synthesised O to H ratio 
should grow with progenitor mass, the significant contribution from 
primordial O will damp the impact of this on the [O i] 6300, 6364 
lines. It is surprising that SN 1987A do es not have the highest value 
since its [O i] emission was shown bv lKozma & FranssorJ ( Il998al) 
to be dominated by emission from the O-rich zones, which should 
be significantly greater in a 20 M© star compared to a low mass star. 
This suggests that either the mass of synthesised O is smaller than 
previously believed, or that SN 1987 A is unusually H-rich. 

The [Ca n] to [O i] ratio is also shown in Figure[8]and suggests 
that either [Ca n] is being overproduced in the model or [O i] un- 
derproduced. A wide range of observed values is seen for this ratio 
but there is no very noticeable trend with SN properties. There may 
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Figure 9. Ratio of flux in [Fen] to [Can] and [Oi]. Tlie dashed black line shows the output from the model. 
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Figure 10. The fraction of the total luminosity emerging in Ha (solid line), 
[Oi] 6300, 6364 A (dotted Une), [Can] 7291, 7323 A (dashed line), and 
[Fe n] 7155, 7172 A (dash-dot line) in the spectral synthesis model. 



be a weaker [Ca n] to [O i] flux ratio for lower velocity SNe such 
as SNe 2008bk and 2003gd seen at -500-550 d compared with 
the model, which could be explained by [Can] emission doing less 
cooling in the H zones at the lower temperatures (of the order of 



lOOOK) seen in low '^Ni mass events (Jerkstrand et al., in prep.), 
but synthesised [O i] still doing most of the cooling in the O zones 
due to fewer competing cooling lines. However, this is a speculative 
explanation and the high- velocity SN 2008cn and 'normal' velocity 
SN 2007aa also show weak [Can] to [Oi] ratios at ~350 d which 
cannot obviously be explained in this way. Further modelling using 
a variety of velocity and '*Ni mass inputs will be used to explore 
this in more detail. 

The other prominent emission line studied here is the [Fe ii] 
7155, 7172 A doublet, which is detailed in |5.1.4| The ratio of [Fen] 
to Hcf model curve is shown in Figure [8] in good agreement with 
the observed range. There is a trend of increasing flux with time 
of [Fen] with respect to Ho- but no obvious trend between SNe 
with different ^*Ni masses and velocities is seen. The ratio of [Fe ii] 
7155, 7172 A doublet to the [Can] 7291, 7323 A hne is plotted in 
Figure |9] again good agreement between the model and observed 
ratios is seen. The [Fen] to [Can] ratio seems to be relatively con- 
stant as a function of time and no obvious trend with ^*Ni mass 
or velocity can be observed. The ratio of the [Fen] 7155, 7172 A 
doublet to the [O i] 6300, 6364 A doublet can be formed using com- 
binations of the primordial and synthesised O and Fe components, 
and is shown in Figure|9] The scatter in the observed values for the 
different SNe is quite large compared to the [Fen] to [Can] ratio 
shown above, perhaps due to varying contributions from synthe- 
sised O to the [O i] lines. 

These comparisons are the first time that nebular phase spec- 
tra have been compared in detail and interpreted in the framework 
of a radiative transfer model. SN1987A has been the cornerstone 
of understanding synthesised material ejected in SNe in relation 
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to a progenitor of known mass, but we now have a sample of 
SNe with progenitor constraints and good quality nebular phase 
spectral coverage. No significant trends with progenitor mass are 
seen, with all the SNe in the sample having very similar flux ra- 
tios for the different elements. This suggests that we are in gen- 
eral, seeing emission from similar regions for the different SNe 
in the sample. The observed line widths and line profiles, along 
with the theoret ical interpretation of the line ratios, suggest like 
for SN 1987 A jKozm a & Fr anssonI Il998al ; iLi & McCra'Jll993j : 
III McCrav & Sunva ev 1993b), that [Can] emission is formed in 
the H zones and the [Fen] emission is produced by a mixture of 
primordial and synthesised Fe. 

However, the most important result discovered here using the 
line profiles and widths and also the theoretical interpretation, is 
that /or our 12 Mq model, the [Oi] 6300, 6364 A emission lines 
are formed from nearly equal contributions from primordial O in 
the H-rich zones and synthesised O in the 0-rich zones. Figure 
[TT] shows a plot of the relative contribution from primordial abun- 
dances as a function of time for the [O i], [Ca ii], and [Fe n] lines. At 
earlier times, the [Can] lines are mainly formed from thermal pro- 
cesses, while at later times, there is also a contribution from fluores- 
cence. Both thermal emission and fluorescence occur mainly in the 
H zone at all times. For the [O i] doublet, the division between syn- 
thesised and primordial material is ~50 per cent each at all times. 
With time, the [Fen] 7155, 7172 A doublet becomes steadily more 
dominated by the H zone, as the Fe/He zone rapidl y cools to tem- 
peratur es where infrared cooling dominates ( Ko zma & FranssonI 
Il998bh . This is markedly different to the results obtained for 
SN 1987 A, which showed that the e mission was being produced 
nearly completely in the O-rich zones jKozma & Franssonlll998al : 
iLi & McCravlll 992VSpvromih o & Pintoll99ll) . The reason for this 
is twofold; the mass of synthesised O is ~ 10 times larger in a 20 Mq 
star compared to a 12 Mq star, and the abundance of O in the en- 
velope of SN 1987A was found to be strongly suppressed by CNO 
mixing. 

5.3 Implications for synthesised O mass estimates 

In the past, the [O i] 6300, 6364 A flux has been used to estimate the 
mass of synthe sised O and hence, the m ain sequence mass of the 
progenitor star. lElmhamdi et al] ( 120033) used an equation that cal- 
culated the O mass in SN 1999em by using a ratio with SN 1987A 
involving its O luminosity in t he [Oi] 6300, 6 364 A lines and its 
^''Ni mass during the tail phase. lUomotol ( [r986ll estimated the min- 
imum mass of O needed to produce the [Oi] emission lines during 
the nebular phase, using an equation that depended on the flux of 
the [Oi] lines, the distance to the SN and the temperature of the 
emitting gas. The O mass of SN 1987A has b een modelled in de- 
tail and found to be in the range 1.2-1.5 Mg t Kozma & FranssonI 
ll998al : IChugailll994i) . 

However, unlike SN 1987A, the IIP SNe studied here result 
from the collapse of significantly less massive stars and we see ev- 
idence that the [O i] emission is made up of contributions from pri- 
mordial and synthesised O, and so the [O i] doublet is not a simple, 
direct indicator of the synthesised O mass. Therefore, the O mass 
estimates in the literature for IIP SNe with lo w progenitor mass es- 
timates should be treated with caution (e.g. Andrews et al. 2011; 
Maguireetal]i2 010; Pozzo et al. 2006; Elmhamdi et al. 2003a; 
Elmhamdil2011 ), since they are based on an analysis of SN 1987 A, 
and the true O masses could be only ~50 per cent of these values. 
However, the O mass can still be estimated as long as the flux con- 
tribution from primordial O in the H zones is subtracted off first. 
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Figure 11. Fraction of [Oi] 6300, 6364 A (dotted line), [Can] 7291, 7323 A 
(dashed line), and [Feii] 7155, 7172 A (dash-dot line) lines being produced 
by primordial abundances in the H zones relative to thek total flux as a 
function of time since explosion. 
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Figure 12. Ratio of the flux in the [O I] 6364 A to the 6300 A lines. The 
black dashed line is the theoretical curve for an O number density of 7x lO'" 
cm"' at 100 d. 



An updated estimate of the synthesised O mass of SN 2004et is de- 
termined in Jerkstrand et al. (in prep.) and can be used in the future 
for comparisons with SNe with low mass progenitor estimates. 



5.4 [O i] 6300, 6364 A line ratio 

The ratio of the 6364 A to the 6300 A line in the [O i] doublet can be 
used to determine whether the O-emitting zones are optically thic k 
or optically thin dSpvromilio & Pintdll99ll ; iLi & McCravHl992h . 
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In Figure [T2I the ratio of tiie 6364 A to the 6300 A hne is shown 
for a sample of eight IIP SNe. In the optically thick limit, the ratio 
of their intensities is 1:1, while in the thin limit the ratio of the 
intensities of the lines should be 1:3. An evolution from optically 
thick toward optically thin is seen for the sample of SNe in Figure 
[T2I The theoretical curve for an O number density of 7 x 10'" cm"' 
at 100 d is also plotted in Figure [T2I 

A contribution of up to ~20 per cent from the Fei 6361 A line 
to the flux in the [Oi] 636 4 A line was found for SN 1987 A by 
iKozma & Franssoril ( Il998ah . This was investigated for this model, 
but at 400 d, the Fe i 6361 A line was found to be a less than 3 per 
cent contaminant and no other lines were found in this wavelength 
region. Instead, we find that opacity due to other lines may influ- 
ence the [O i] line ratio. The escape probability for photons, with 
respect to absorption in other lines is dependent on wavelength. 
For the two [O i] 6300, 6364 A lines there is a ~60 per cent and 
~75 per cent chance of escaping, respectively. This will increase 
the 6364/6300 ratio, since the 6364 A photons will escape some- 
what more easily (by a factor of 1.25 in this case). 

There is a significant difference in the rate at which the O- 
emitting regions become optically thin for different SNe. SN 2004et 
shows a steep decline towards the optically thin regime, which 
might suggest that the density of its ejecta is low, which can be 
explained by the high HWHM [O i] velocities seen for SN 2004et 
of ~1200 kms~' during the nebular phases studied. In a similar 
way, the low velocity SN 2008bk had lower HWHM [O i] veloc- 
ities (~350 kms"') and therefore remained denser for longer and 
hence took longer for its ejecta to become optically thin. However, 
SN 1987 A like SN 2004et, had high [Oi] velocities (~1 100 km s"') 
and therefore, should have become optically thin at a similar rate to 
SN 2004et, but this evolution is not seen. This could be due to the 
large contribution from the synthesised [O i] in the core, where the 
densities are higher. 

Indeed, this ratio can be used to probe the density struc- 
ture of the O-emitt ing regions as was done for SN 1987 A by 
iLi & McCravl iwj) and'Spvromilio & Pintd h99lh . However, un- 
like SN I987A, it has been shown for lower mass stars, that the [O i] 
6300, 6364 A doublet is formed equally by contributions from the 
envelope and 0-rich zones and so any density estimates must use a 
two component model to fit the contributions from both the H-rich 
and O-rich zones. 



6 LINE IDENTIFICATIONS OF SN 2008BK 

A late-time nebular spectrum (-1-547 d) of SN 2008bk is shown in 
Figure [T3] along with a model spectrum outputted from the spectral 
synthesis code. The identifications do not distinguish between for- 
bidden and allowed lines, instead it identifies which element con- 
tributes most to which features. This SN has very low line velocities 
that gives narrow emission lines that are well resolved. In higher 
velocity events, such as SN 2004et, many lines are blended not due 
to the resolution limit of the instruments used to obtain the spec- 
tra but because their profiles intrinsically overlap in velocity space. 
Therefore, low velocity SNe give us a great opportunity to identify 
individual line emissions that would be blended in higher velocity 
events. The prominent emission lines in this spectrum were iden- 
tified using the same spectral synthesis code as discussed earlier, 
again using a 12 M© model, but this time choosing a core veloc- 
ity of 600 km s"' to match the velocities seen in SN 2008bk. Adi- 
abatic is strong relative to line cooling in the envelope, as the low 
^^Ni mass implies a low electron fraction and low collision frequen- 



cies. This makes the steady-state temperature calculation unreliable 
and therefore, this model has not been used to calculate flux ratios. 
However, the line identifications should not be affected and can be 
used to identify the spectral features in low velocity spectra such as 
those of SN 2008bk. 

Most of the lines in the complex region from 4000-5500 A 
can be identified with [Fei] and [Feii] lines, created by scattering 
and fluorescence in the H envelope. Above 5500 A the core zones 
start making contributions to the spectrum. The 5900 A line is Nai 
5896 A which scatters emission from the Nai 5890 A line, as well 
as from the Hei 5876 A line, with about a third of the photons 
emitted by synthesised Na in the 0/Ne/Mg zone. The Fe i lines at 
7900, 8030, 8060, 8200, 8300, and 8350 A are all from Fe in the 
core (both from the Fe/He clumps and the 0/Ne/Mg clumps). In the 
near-infrared, we can identify the O i 8446 A line, which is usually 
hard to see due to blending with the Ca n triplet in higher- velocity 
SNe. Another line which usually also suffers from this blending is 
[Ci] 8727 A which is clearly seen here. In the model, this line, 
along with the Ci lines at 9000-9100 A and 9400 A are created 
by C in the He envelope, verifying that significant dredge-up has 
occurred in this star. 



7 CORRELATION OF VELOCITY WITH ^'^Ni MASS 

There are currently a number of methods by which the ^*Ni mass 
ejected in a IIP SN can be estimated. Most methods require knowl- 
edge of the bolometric light curve of the SN during the early ra- 
dioactive tail phase. During this phase, the bolometric light curve 
is powered nearly completely by the radioactive decay of ^^*Co to 
'*Fe and hence the luminosity at this stage gives a good indication 
of the synthesised ^*Ni mass that is ejected in the explosion. The 
link between the light curve and ^*Ni mass can be made in three dif- 
ferent ways: a comparison with the luminosity of SN 1987A, which 
ha s a well constrained '*Ni mass, t he stee pness parameter method 
of lElmham di. Chugai & D anzigej ( l2003bl) , and the tail luminosity 
method of iHamuvi t2003) . However, the use of these methods is 
limited to SNe with good bolometric coverage during the radioac- 
tive tail phase, as well as reliable distance and extinction es timates. 

For S Ne with limited multi- band coverage, Berst en & Hamuvl 

( I2OO9I) and lMaguire et al.l j2010l) described bolometric corrections 
to convert from a magnitude obtained in a single band to a bolo- 
metric luminosity. The former paper estimated the correction as a 
function of colour, while the latter work as a function of time since 
explosion. These corrections are useful for sparsely covered events, 
but still require an accurate estimate of the extinction and distance 
to the S N. Us ing late-tim e spectra (200-400 d po s t-explos ion), 
IChugail ( Il990h and Elmh amdi. Chugai & DanzigeJ j2003bh de- 
scribed how the luminosity of the Ho- emission line is proportional 
to the ejected ^^Ni mass. This method is useful when no photom- 
etry is available, but again the distance and extinction to the SN 
must be known. Therefore, independent methods that do not rely 
on knowledge of the ab solute magnitude of the SN would be a 
useful addition. Hamuvl ( l2003h showed a correlation between the 
'*Ni mass and the photospheric velocity at mid-plateau. However, 
the photospheric velocity decreases rapidly during the photospheric 
phase and therefore, requires a well-sampled light curve to provide 
a secure constraint on the explosion epoch to obtain the velocity 
mid-plateau. In this section, we discuss a method for estimating the 
ejected '*Ni mass from the FWHM of the emission lines observed 
in nebular phase spectra of IIP SNe, where the velocity is found 
to relatively constant in the range ~350-550 d post explosion and 
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Figure 13. Line identification for tlie low-velocity, SN 2008bk at +547 d post-explosion comparing the spectral synthesis model (blue) to the data (red). We 
have not distinguished between forbidden and allowed lines, instead we show which ion contributes most to which feature. 
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Figure 14. The FWHMcon of the Ha emission line corrected for the in- 
strumental resolution for SNe 1987A, 1997D, 1999em, 2003gd, 2004et, 
2008bk, and 2009N as a function of time since explosion. 

hence does not require accurate constraints on the explosion epoch. 

The sample of SN spectra used in this work is detailed in Sec- 
tion[2] The sample is a combination of IIP SNe found in the litera- 
ture with good photometric and spectral coverage during the nebu- 
lar phase, along with new previously unpublished observations. The 
photometric coverage was necessary so that the '^Ni mass could be 



estimated independently from the spectra and an empirical relation 
between the ejected ^^Ni mass and the line widths could be deter- 
mined. We have added one more IIP SNe with nebular phase obser- 
vations and reliable '*Ni mass estimates to the sample, SN I997D, 
which was a low-l uminosity, low velocit y SN with an ejected '*Ni 
mass of 0.008 Mq jPastorello et al.l2004h . SN 1997D is not studied 
in detail here because of the lack of informat ion on its progenitor 
star. The spectra of SN 1997D are taken from lBenetti et allllOOlh 
with a spectral resolution of 13 A. 

The FWHM of the most prominent feature in the spectra, the 
Ho- 6563 A emission line was measured by fitting a Gaussian to the 
emission profile. This line is very strong in all the spectra and is 
unblended with other lines so we can be confident that the FWHM 
can be measured to good accuracy. As in Section |4~2l the FWHM 
of the line profiles were corrected for the spectral resolution of 
the instruments used to obtain each spectrum. Figure[74lshows the 
FWHMcorr of the Ha emission line as a function of time for 8 IIP 
SNe that also have ^^Ni masses estimated from their bolometric 
light curves. From Figure [141 it can be seen clearly that for any 
given SN, the FWHMcon is relatively constant at these epochs and 
so observations at a coeval epoch are unnecessary for a compar- 
ison to be made. A mean value of the FWHMcon- was calculated 
for each SN from the available data at epochs in the range ~350- 
550 d. These values can then be compared to the ^*Ni mass esti- 
mated from the luminosity of the nebular phase light curve. The 
FWHMcon value used for each SN along with its ^^Ni mass esti- 
mate from the literature are given in Table |5] 

Figure[T5]shows the ejected ^*Ni mass calculated from the tail 
phase luminosity of the SNe against the FWHMcon of the Ha emis- 
sion line. It is immediately apparent that a correlation exists be- 
tween these parameters, with small FWHMcon (low velocity) events 
such as SNe 1997D, 2008bk, and 2009N displaying correspond- 
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Figure 15. The ^''Ni mass obtained from the tail phase of the bolometric 
light curve against the FWHM of the Ha emission line for the SNe listed in 
Tabled 



these SNe. The uncertainties in these estimates are relatively large, 
but at the very least this method can be used to determine between 
low, intermediate and high ejected ^*Ni masses. A larger sample is 
needed to study this method of ""^Ni mass determination in more 
detail and in particular it would be of interest to include some IIL 
SNe to see if the relation also holds for these events. 

A correlation bet ween the ^^Ni m ass and explosion energy for 
IIP SNe was shown bv lHamuvl ( |2003|) for a sample of 16 IIP SNe, 
with higher ^^Ni mass events having greater explosion energies. 
iNadvozMi] ( |2003|) also describe a correlation between the explo- 
sion energy and the ^^Ni mass of IIP SNe, which was suggested to 
have a complex nature that also depends on the progentior mass. 
One can assume that the velocity (or equivalently the FWHM) of 
the line profiles will scale with the explosion energy (which is a 
combination of internal and kinetic energy), since the more en- 
ergy produced in the explosion, the more kinetic energy that will 
be given _to t he e jecta (at least 90 per cent of the explosion en- 
ergy tArnetllil99d) ). This is the first time a correlation between 
the FWHM of the nebular phase line profiles and the ^^Ni mass 
has been identified, and assuming similar ejecta masses for the SN 
sample, can be interpreted as a correlation between the explosion 
energy and the '^Ni mass. 



Table 5. Data used for Figure [Tll to determine the relationship between the 
ejected ^^Ni mass and the FWHM of Ho- corrected for the spectral resolu- 
tion. 



SN 


FWHMhq, (A) 


5*Ni mass (Mo) 


Reference 


1997D 


20 ±4 


0.008 ± 0.006 


Benettietal. (2001) 


2008bk 


21 ±2 


0.004 ± 0.003 


Pignata et al. (in prep.) 


2009N 


28 ±3 


0.009 ± 0.004 


Takats et al. (in prep.) 


2003gd 


44 ±5 


0.016 ±0.01 


Hendrv et al. (2005) 


I999em 


51 ±2 


0.036 ± 0.005 


Utrobin (2007) 


2004et 


60 ±2 


0.064 ±0.04 


Maeuire etal. (2010) 


1987A 


70 ±4 


0.069 ± 0.0012 


Bouchetetal. (1991) 



ingly smaller Ni masses. In a similar way, the higher velocity 
SNe of the sample, such as SNe 1987A and 2004et, have higher 
^''Ni masses. The relationship between the mass of ejected ^*Ni and 
the FWHMcorr of the Ha emission line can be fitted with an equa- 
tion of the form, 

M(^^^Ni) = A X 10'^ FWHM„„ 

where B = 0.0233 ± 0.0041 and A = 1.81!^°^ x 10^^ The re- 
duced value for this fit is 0.73. This equation can used for SNe 
for which limited photometric data are available or where there are 
large uncertainties on the distance and extinction estimates. 

As an example we apply this relation to two IIP SNe, with- 
out good late-time light curves. SN 2008cn is a IIP SN that was 
observed spectroscopically at the VLT at an epoch of +358 d post- 
explosion. No photometric data has yet been published on this ob- 
ject, but from our spectrum, using a FWHMemr of Ha of 79 A, we 
have estimated the ^*Ni mass to be 0.13!° ^" Mq. SN 2007aa is an- 
other IIP SN that was observed at the VLT at a nebular phase epoch 
of -1-376 d post-explosion and found to have for Ha a FWHM^on 
of 52 A which gives a '''Ni mass of 0.03;^q[5, Mq. Hopefully these 
'*Ni mass estimates will be confirmed by future data releases for 



8 SUMMARY AND CONCLUSIONS 

In this paper we have investigated the nebular phase spectra of a 
sample of IIP SNe and compared their observational analysis to a 
detailed radiative transfer model, to try to better understand the for- 
mation of their emission lines and the variation within the sample. 
We have focused on IIP SNe with either progenitor detections or 
restrictive limits on their progenitor masses. Despite the large vari- 
ation in '*Ni masses and expansion velocities, the characteristics of 
the emission lines in the nebular spectra of IIP SNe are found to be 
relatively homogeneous and independent of progenitor masses. For 
each SN, the shapes of the different emission line profiles are very 
similar and suggest a similar spatial distribution of the different ele- 
ments. This is also seen when the velocities of the prominent emis- 
sion lines are measured using their HWHM, with only Ha showing 
a noticeably higher velocity, which suggests a larger contribution 
from farther out in the ejecta. 

For SN I987A, as has been previously shown, we have shown 
that the Ha velocity was significantly larger than the [O i] veloc- 
ity, while for the other SNe in our sample the differences between 
H and [O i] were not as pronounced. This suggests that the zones 
responsible for [Oi] emission are different in the rest of the SN 
sample compared to those in SN 1987 A. From analysis of the line 
profiles and widths, it appears that all IIP SNe (except SN I987A) 
have similar spatial distributions of [O i] which is closer in veloc- 
ity to Ha than for SN I987A. This implies that the [Oi] flux could 
be predominantly emitted by primordial O in the H zones and may 
not be coming from the synthesised O from the inner regions of the 
ejecta. The line profiles are also shown to be intrinsically peaked in 
shape, even after taking into account the convolution with the in- 
strumental resolution. This shows that mixing must have occurred 
for all the elements to be producing emission from both the centre 
and the outer regions of the ejecta. 

These observational results are corroborated by the output of 
the radiative transfer model of a 12 Mq star, where a significant 
proportion of the flux of the prominent emission lines is found to 
originate in the H-rich zones. For SN I987A, the [Ca ii] and several 
[Fe ii] emission lines were also found to be predominantly formed 
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in the H zones, but here we present the new result that for the 12 Mq 
progenitor model used here, the [O i] flux is produced in approxi- 
mately equal amounts by O in the H zones and O in the inner O- 
rich zones, and this explains the much smaller difference between 
the Ha and [Oi] velocities seen for our sample of SNe. This has 
major implications for current methods of calculating the synthe- 
sised O mass in IIP SNe, which assume that the [O i] 6300, 6364 A 
line emission is produced solely by synthesised O. A calculation of 
the synthesised O mass of SN 2004et, that takes into account these 
new results will be presented in Jerkstrand et al. (in prep.). 

The flux ratios of the emission lines from the 12 M© model 
were found in general to agree quite well with the observed flux 
ratios. No trends were found between progenitor mass and the ob- 
served flux ratios, in agreement with the result that the flux in these 
lines has a significant contribution from primordial abundances. We 
also investigated possible correlations between the spectral prop- 
erties and other observational properties of the SNe and found a 
correlation between the ejected ^^^Ni mass and FWHM of the Ha 
emission line. An empirical relation was obtained, which can be 
used to estimate the ^*Ni mass of SNe with limited bolometric cov- 
erage and/or poor constraints on their extinction and distance. The 
benefit of this method is that it requires only the width of the Ha 
emission line during the nebular phase to estimate the '*Ni mass. 
Even more importantly, this provides a strong constraint on explo- 
sion models, particularly the amount of fall-back material in them. 
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